Analyzing 77' photoproduction data on the proton at energies of 1.5—2.3 GeV 
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The recent high-precision data for the reaction 7p pr]' at photon energies in the range 1.5- 
2.3 GeV obtained by the CLAS collaboration at the Jefferson Laboratory have been analyzed within 
an extended version of the photoproduction model developed previously by the authors based on a 
relativistic meson-exchange model of hadronic interactions [Phys. Rev. C 69, 065212 (2004)]. The 
T)' photoproduction can be described quite well over the entire energy range of available data by 
considering Sii, Pii, Pia, and D13 resonances, in addition to the t-channel mesonic currents. The 
observed angular distribution is due to the interference between the f-channel and the nucleon s- 
and u-channel resonance contributions. The j — 3/2 resonances are required to reproduce some of 
the details of the measured angular distribution. For the resonances considered, our analysis yields 
mass values compatible with those advocated by the Particle Data Group. We emphasize, however, 
that cross-section data alone are unable to pin down the resonance parameters and it is shown that 
the beam and/or target asymmetries impose more stringent constraints on these parameter values. 
It is found that the nucleonic current is relatively small and that the NNrj' couphng constant is not 
expected to be much larger than 2. 

PACS numbers: 25.20.Lj, 13.60.Le, 14.20.Gk 



I. INTRODUCTION 

One of the primary interests in investigating the rj' 
photoproduction reaction is that it may be suited to ex- 
tract information on nucleon resonances, N*, in the less 
explored higher N* mass region. Current knowledge of 
most of the nucleon resonances is mainly due to the study 
of TriV scattering and/or pion photoproduction off the nu- 
cleon. Since the rj' meson is much heavier than a pion, 
rj' meson-production processes near threshold necessar- 
ily sample a much higher resonance-mass region than 
the corresponding pion production processes. They are 
well-suited, therefore, for investigating high-mass reso- 
nances in low partial-wave states. Furthermore, reaction 
processes such as 77' photoproduction provide opportuni- 
ties to study those resonances that couple only weakly 
to pions, in particular, those referred to as "missing res- 
onances" , which are predicted by quark models, but not 
found in more traditional pion-production reactions [1] . 

Another special interest in 77' photoproduction is the 
possibility to impose a more stringent constraint on its 
yet poorly known coupling strength to the nucleon. This 
has attracted much attention in connection with the so- 
called "nucleon-spin crisis" in polarized deep inelastic 
lepton scattering [2]. In the zero-squared- momentum 
limit, the NNrf coupling constant gNNrj'iq'^ = 0) is re- 
lated to the flavor-singlet axial charge Ga through the 
flavor singlet Goldberger-Treiman relation [3] (see also 
Refs. [4-6]) 

2mN Ga{0) = F givNrt'iO) + TTTT ml gNNoiO) , (1) 
where is a renormalization-group invariant decay con- 



stant defined in Ref. [3],^ Np is the number of flavors, 
and mjv and m^/ are the nucleon and rj' masses, respec- 
tively; gNNG describes the coupling of the nucleon to the 
gluons arising from contributions violating the Okubo- 
Zweig-Iizuka (OZI) rule [7]. The EMC collaboration [2] 
has measured an unexpectedly small value of Ga(0) ~ 
0.20 ±0.35; a more recent analysis of the SMC collabora- 
tion [8] yields a comparable value of Ga{Q) ~ 0.16 ±0.10. 
The first term on the right-hand side of the above equa- 
tion corresponds to the quark contribution to the "spin" 
of the proton, and the second term to the gluon contri- 
bution [5, 9]. Therefore, once gNNr)'{0) is known, Eq. (1) 
may be used to extract the coupling gNNciO)- Unfortu- 
nately, however, there is no direct experimental measure- 
ment of gNNri'{0) so far. Reaction processes where the 77' 
meson is produced directly off a nucleon may thus offer 
a unique opportunity to extract this coupling constant. 
Here it should be emphasized that, as has been pointed 
out in Ref. [10], hadronic model calculations such as the 
present one cannot determine the NNrj' coupling con- 
stant in a model-independent way. At best, we get an es- 
timate for the range of its value at the on-shell kinematic 
point, i.e., at = rn^,. Assuming the usual behavior of 
hadronic form factors for off-shell mesons which generally 
decrease for < m^, we expect then that an eventually 
small upper limit of gNN-q'iq'^ = ™^') would lead to an 
even smaller value of gNNrj' (0), which is needed in Eq. (1) 
to extract gNNciO)- 



1 In the OZI limit, F 
decay constant. 



/2NpF,r, where stands for the pion 
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The major purpose of the present work is to perform an 
analysis of the 7p prj' reaction within an extended ver- 
sion of the relativistic meson-exchange model of hadronic 
interactions as reported in Ref. [10]. This analysis is mo- 
tivated by the new high-precision cross-section data ob- 
tained by the CLAS collaboration [11] at the Jefferson 
Laboratory (JLab). The new data supersede the previ- 
ous SAPHIR data [12] analyzed in Ref. [10] both in ab- 
solute normalization and angular shape. Also, the new 
CLAS data are much more accurate and, as such, may 
reveal features that were not seen in the analysis of the 
SAPHIR data. 

The present paper is organized as follows. In Sec. II 
the extension of our model [10] for jp prj' is given. 
The results of the corresponding model calculations are 
presented in Sec. III. Section IV contains a summary with 
our conclusions. Some technical details of the present 
model are given in the Appendix. 

II. FORMALISM 

The dynamical content of the present rf photoproduc- 
tion calculation is summarized by the graphs of Fig. 1 
where we employ form factors at the vertices to account 
for the hadronic structure. The gauge invariance of this 
production current is ensured by a phenomenological con- 
tact current, according to the prescription of Refs. [IS- 
IS]. This contact term provides a rough phenomenologi- 
cal description of the final-state interaction which is not 
treated explicitly here. The basic details of the present 
approach are the same as in our previous paper [10] and 
we will not repeat them here. There are, however, a few 
improvements and those will be discussed here. 




FIG. 1: Diagrams contributing to 7p — » rj'p- Time proceeds 
from right to left. The intermediate baryon states are denoted 
l\l for the nucleon, and R for the Sn, Pn, P13, and D\3 reso- 
nances. The intermediate mesons in the t-channcl arc p and 
u). The external legs are labeled by the four-momenta of the 
respective particles and the labels s, u, and t of the hadronic 
vertices correspond to the off-shell Mandelstam variables of 
the respective intermediate particles. The three diagrams in 
the lower part of the diagram are transverse individually; the 
three diagrams in the upper part are made gauge-invariant by 
an appropriate choice for the contact current depicted in the 
top-right diagram. The nucleonic current (nuc) referred to in 
the text corresponds to the top line of diagrams; the meson- 
exchange current (mec) and resonance current contributions 
correspond, respectively, to the leftmost diagram and the two 
diagrams on the right of the bottom line of diagrams. 

transition current reads 

_ ^92m^ / q^n*Pq y^±)n) F'^^ + H. c. , 

(4) 

where F'"' = d^A'' — d" A^^ is the electromagnetic field- 
strength tensor (with A^ being the vector potential). 



A. Spin-3/2 resonances 

The present fits also require the inclusion of spin-3/2 

resonances, denoted generically by N* . The Lagrangian 
for the hadronic NN*rj' interaction is given by 

£W = ^^^m^*/.0^^(^)r(±)jva- V + H. c. , (2) 

where iV*'^, N , and 77' are the resonance, nucleon, and 
meson fields, respectively, and 

r(+) = 75 and F^-) = 1 (3) 

pertain to positive- and negative-parity resonances, re- 
spectively. For the couphng tensor, = ^f^^ ~ {z + 
i)7p7,y, we take z = for the off-shell parameter 
for simplicity.^ The Lagrangian for the electromagnetic 



B. Energy-dependent resonance widths 

For the present application, we have adapted our 
formalism to accommodate energy-dependent resonance 
widths with the appropriate threshold behavior. 

For a spin- 1/2 resonance propagator, we use the ansatz 

a / -V ^ 1 ^ i) + mR 

(5) 

where mR is the mass of the resonance with four- 
momentum p. r is the width function whose functional 
behavior will be given below. 

For spin-3/2, the resonant propagator reads in a 



^ We have also explored how the fits changed upon varying the off- 
shell parameter z in Eqs. (2) and (4), since spin-3/2 resonances 



also play a relevant role in reproducing the data quantitatively, as 
discussed in Sec. III. However, we didn't observe any significant 
changes due to this parameter and we feel justified, therefore, to 
keep this parameter at z = — | for simplicity. 
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schematic matrix notation 

'53/2 (P) 



(j/)-mR)g - i—T 



A . 



(6) 



All indices are suppressed here, i.e., g is the metric ten- 
sor and A is the Rarita-Schwingcr tensor written in full 
detail as 



1 



3m\ 



(7) 

where /?, a, and e enumerate the four indices of the 7- 
matrix components (summation over e is implied). The 
inversion in (6) is to be understood on the full 16- 
dimcnsional space of the four Lorcntz indices and the 
four components of the gamma matrices. The motiva- 
tion for the ansatz (6) and the technical details how to 
perform this inversion is given in the Appendix. 

In both cases, we write the width F as a function of 
W = -v/s according to 



T{W) = Tr 



N-, 



(8) 



where the sums over i and j respectively account for de- 
cays of the resonance into N two- or three-hadron chan- 
nels and into radiative decay channels. The total 
static resonance width is denoted by Tr and the numer- 
ical factors l3i and (with < /3j,7j < 1) describe the 
branching ratios into the various decay channels, i.e., 



JV 



(9) 



3=1 



Similar to Refs. [16-19], we parameterize the width func- 
tions Ti and T^^ (which arc both normalized to unity 
at W = tur) to provide the correct respective threshold 
behaviors. 

For the decay of the resonance into two hadronic frag- 
ments with masses mn and mj2, the hadronic width func- 
tions Vi are taken as 



ti{W) = 



_cu_ 



Di{W) (10) 



for W > m,ii + rrii2, and zero otherwise. L denotes the 
partial wave in which the resonance is found and the mo- 
mentum qj is the magnitude of the center-of-momentum 
three-momentum of the two fragments, i.e.. 



C^^{W) 



^/[W^ - {ma + rrii^y] [W^ - (m^ - mia)^] 



2W 



(11) 

and qiR — qilniR). For the decay of the resonance into 
one baryon and two mesons (for example, Ntttt), we use 



ti{W) 



(liR 



2L+4 



Af+q^ 



iR 



Af + qf 



L+2 



Di{W) , (12) 



where in (11) needs to be replaced by the sum of the 
two meson masses for this case, and mn is the baryon 
mass. In principle, the factor 



Di{W) 



\W J 



with 



> 



(13) 



allows for a modification of the asymptotic behavior of 
Ti{W), however, we use rij = throughout for simplic- 
ity. The parameter Aj is an inverse range parameter; 
since we found very little sensitivity to varying this pa- 
rameter (within reasonable ranges), we kept it fixed at 
Aj = 1 fm~^ for all channels. 

The width function F^. for the decay into a hadron 
with mass rrij and a photon with three-momentum kj is 
taken as 

where 



2L+2 



L+l 



(14) 



2W 



(15) 



for W > rrij, and zero otherwise, and kjR = kj{mR). As 
in the hadronic case, the asymptotic damping function is 
given by 



with r,- > 



(16) 



Again, for simplicity, we employ rj = throughout. In 

practice, for the present case, the photon decay chan- 
nels are negligibly small and play no role for the total 
width. The corresponding branching ratio jnj for the 
Nj channel is only needed to extract the value of the 
Nr]' branching ratio Pnt]' (see below). 



III. RESULTS AND DISCUSSION 

Before we discuss the details of our results, some gen- 
eral remarks are in order. The basic strategy of our 

model approach is to start with the nucleon plus meson- 
exchange currents and add the resonances one by one as 
needed in the fitting procedure until one achieves a rea- 
sonable fit of the new 77' photoproduction data obtained 
by the CLAS collaboration [11]. We allow for both spin- 
1/2 and -3/2 resonances in our model. Our quantitative 
criterion for a reasonable fit was to discard all fits with 
a per data point oi /N > 1.3, which is supported 
by the fact that fits with 'y^ jN much larger than 1.3 are 
noticeably of inferior fit quality even for the naked eye. 
Under this criterion, we found that one needs at least 
four resonances in order to obtain a reasonable fit in the 
present approach. We find, in particular, that, in addi- 
tion to spin- 1/2 resonances, spin-3/2 resonances are nec- 
essary to achieve acceptable fits. In this respect, we em- 
phasize that the SAPHIR data [12] analyzed in Ref. [10] 
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have rather large error bars. While not entirely incom- 
patible with the new high-precision CLAS data [11], they 
clearly are less constraining than the CLAS data, which 
may explain why there was no need for spin-3/2 reso- 
nances in our previous work. 

As we have pointed out in Ref. [10], the cross-section 
data alone are unable to pin down the model parame- 
ters and, therefore, one finds different sets of parameters 
which fit the data equally well. Note that this is not 
due to the uncertainties in the data, but simply because, 
intrinsically, the cross sections do not impose enough 
stringent constraints on the fit. In particular, for each 
resonance, the resulting fitted mass value depends to a 
certain extent on its starting value in the fitting proce- 
dure. The starting (resonance) mass values we consider 
here generally are around those advocated by the Particle 
Data Group (PDG) [20]. 

In the present work, in the case of those resonances 
that can be identified with known PDG resonances, we 
have taken into account only the corresponding domi- 
nant branching ratios /?; from the PDG for hadronic de- 
cays when this information is available (and we ignored 
the fact that some of the quoted branching ratios are 
subject to large uncertainties). If no information is avail- 
able, we consider only the Nit partial decay, with the 
corresponding branching ratio (iisi-n as a free fit parame- 
ter. Apart from these branching ratios, we also consider 
the Nr]' branching ratio l3Nr)' which is calculated from 
the product of the coupling constants gNN'-q'SNN'-j in 
conjunction with the assumed branching ratio ^M-y for 
the radiative decay. In the following tables, therefore, 
Pnyi' is not an independent fit parameter, but rather a 
parameter extracted from the fitted values of the product 

9NN*ri'gNN'i- 

One might expect that the way in which the energy de- 
pendence is implemented in the resonance width in the 
present work [cf. Eqs. (12)-(14)] may introduce a consid- 
erable uncertainty in the final results. However, we find 
that the cross sections are not very sensitive to our as- 
sumptions in this respect. In fact, we also re-ran some of 
the cross-section fits discussed below using step functions 
for the widths that switch on the full partial widths at the 
corresponding thresholds without any smooth energy de- 
pendence and we found that the parameter sets obtained 
in this way were fairly close to the ones reported here. 
For spin observables, however, this insensitivity does not 
hold true. In particular, the beam and target asymme- 
tries are rather sensitive to how the energy dependence of 
the width is treated and one must be careful then when 
confronting model predictions with the data when the 
latter should become available. 

We now turn to the discussion of the details of our anal- 
ysis. We emphasize that the results shown here do not 
necessarily have the lowest . Rather, they are sample 
fit results that illustrate the different dynamical features 
one may obtain considering only the currently available 
data in the analysis within the fit-quality criteria men- 
tioned above. 
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FIG. 2: (Color online) Differential cross section for 7p — > pr]' 
according to the mechanisms shown in Fig. 1 as a function 
of the r;' emission angle in the center-of-momentum frame 
of the system. As indicated in the legend, the curves corre- 
spond to the fit results of Tables I (dash-double-dotted), II 
(dotted lines). III (dashed), IV (solid), and V (dash-dotted). 
The numbers (T7, W) in parentheses are the incident photon 
energy T-y and the corresponding s-channel energy W — y/s, 
respectively, in GeV. The data are from Ref. [11]. 



A. Differential cross sections 

The details of the fits presented here are given in Ta- 
bles I-V and the corresponding Figs. 3-7. For the pur- 
pose of easy comparison, Fig. 2 provides an overview of 
all those results. All five fits were obtained using the 
energy-dependent width functions of Sec. II B. They 
all have comparable overall from each other and de- 
scribe the data quite well. We see that most of the dif- 
ferences among them are at forward and backward an- 
gles where there are no data. Therefore, measurements 
of the cross sections at more forward and backward an- 
gles than presently available would disentangle some of 
the results in Fig. 2. Despite the fact that the over- 
all quality of the fits is comparable to each other, the 
resulting parameter values are quite different. In partic- 
ular, the fit set in Table I contains the minimum num- 
ber of resonances (four) required to meet the present fit- 
quality criterium mentioned in the beginning of this sec- 
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TABLE L Model parameters fitted to the 7p — + rj'p. (See text, 
and also Ref. [10], for explanations of parameters.) Values in 
boldface are not fitted. The branching ratios 7JV7 are assump- 
tions made to extract /3jv^' (for the total width, however, the 
7iV7 values are too small to be relevant). The starting values 
for fitting all resonance masses were chosen here within the 
energy range covered by the data set. "x^ jN — 1.19. 



Nucleonic current: 




dNNri' 


0.43 


X 


0.0 


An (JVleV) 


1200 


Mesonic current: 






1.25 




0.44 


A„ (MeV) 


1275 


N* — Su current: 




mjv (MeV) 


1958 


<?JVJV*7 gNN'-qi 


0.25 


X 


1.00 


A / TV T t 7'\ 

Ajy. (MeV) 


1200 


Vn* (MeV) 


139 




0.002 




0.50 


pNn' 


0.50 


N* — Pii current: 




JTljv* (MeVj 




QNN'-f QNN'ri' 


0.80 


X 


1.00 


An* (MeV) 


1200 


Vn* (MeV) 


136 


7iV7 


0.002 




0.36 





U.D4 


N* = P13 current: 




jTijv (MeV) 


1885 


giNN'j gNN'T]' 


0.01 




0.17 


A /'l\/T^a\/^ 

iViv* t^ivie V ) 




rjv- (MeV) 


59 


/3jV7r 


0.6 




0.4 


N* — D13 current: 




TUN' (MeV) 


1823 


giNN*-t gNN'ri' 


0.47 


g2NN*i gNN'r]' 


-0.65 


rjv- (MeV) 


450 


7iV7 


0.002 




1.00 



tion. In contrast to the analysis of the SAPHIR data 
[10], the inclusion of the spin-3/2 resonaces is impor- 
tant in order to reproduce the data quantitatively. As 
mentioned before, although we cannot identify the reso- 
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FIG. 3: (Color online) Differential cross sections and the dy- 
namical content of the present model corresponding to the 
fit result of Table I. The dash-dotted curves correspond to 
the mesonic current contribution; the dashed curves to the 
511 resonance current and the dotted curves to the P\\ reso- 
nance. The dot-double-dashed curves correspond to the P13 
resonance current while the dash-double-dotted curves show 
the D13 resonance contribution. The solid curves correspond 
to the total current. The nucleonic current contribution (long- 
dashed curves) is very small and cannot be seen on the present 
scale. 



nances uniquely in the present analysis. Table I reveals 
that one of the resulting resonances, Pii(2104), is consis- 
tent with that quoted by the PDG [20] as one-star reso- 
nance. In the fit set of Table II, we have included an ad- 
ditional D13 resonance. Here, all the resonances but one 
are above the 77' production threshold energy and that 
two of the resulting spin-3/2 resonances, ^13(1900) and 
£'13(2084), are consistent with those seen and quoted by 
the PDG [20] as two-star resonances. Also, in this par- 
ticular set of parameters, the resulting NNrj' coupling 
constant is very small. The fit set of Table III includes 
three 5*1 1 and three Pn resonances, instead of one each as 
in Table II, keeping the number of spin-3/2 resonances 
unchanged compared to the fit set of Table II. Here, 
two of the Sii, one of the Pn and one of the D13 reso- 
nances end up well below the production threshold, while 
one Pii resonance mass is close to 2.4 GeV. With the 
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TABLE 11: Same as Table. I. An additional Dia resonance 
was included here to see whether this would improve the fit 
quality. xV^ = 1-04. 



Nucleonic current: 






SnNt}' 


0.25 




A 


U.U 




A /TV yf T 

An (MeV) 


1200 




Mesonic current: 






dr}' PI 


1.25 






0.44 




A„ (MeV) 


1308 




A^* = Sii current: 






mjv. (MeV) 


1925 




gNN'-i gNN'ri' 


0.08 




A 


0.08 




Ajv* (MeV) 


1200 




Tjv. (MeV) 


40 




7iV7 


0.002 




I3ntt 


0.56 






0.44 




N* — Pii current: 






mjv* (MeV) 


1991 




gNN'-f QNN'ri' 


-1.69 




A 


0.09 




Ajv* (MeV ) 


1200 




Tjv. (MeV) 


158 




7JV7 


0.002 






0.42 






0.58 




N* — Pi3 current: 






mjv* (MeV) 


1907 




giNN'-f gNN'ri' 


-0.06 




g2NN'f gNN'ri' 


-0.09 




Ajv* (MeV) 


1200 




1 jv (MeV) 


123 




7JV7 


0.002 




/3jV7r 


0.60 






0.4 




Pnti' 


0.00 




N* — Di3 current: 






771 A^* (^iviev j 




zUo4 


giNN'j gNN'ri' 


-1.17 


-0.21 


g2NN'-t gNN'ri' 


0.53 


0.19 


An' (MeV) 


1200 


1200 


rjv* (MeV) 


55 


108 


7JV7 


0.002 


0.002 


Pntt 


1.00 


0.54 


Pnt)' 


0.00 


0.46 
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FIG. 4: (Color online) Same as in Fig. 3 for the fit result 
of Table II. The nucleonic current contribution (long-dashed 
curves) is very small and cannot be seen on the present scale. 



PDG [20] as four-star [511(1535), S'ii(1650)], three-star 
[Pii(1710),Di3(1720)], two-star [Pi3(1900), i^i3(2080)], 
and one-star [S'ii(2090), Pii(2100)] resonances. Here, 
the NNt]' couphng constant is gNNri' ~ 1-3. In the 
fit result of Table IV, we have omitted the P13 reso- 
nance and considered two 5'ii, three Pn and two I?i3 
resonances. Again, three of the resulting resonances, 
S'ii(1538), Pii(1710), and i:)i3(2090), are consistent with 
known resonances. The NNrj' coupling constant is prac- 
tically zero, in line with the small value obtained for the 
fit result of Table II. We have also considered all the 
known spin-1/2 and -3/2 resonances [20] (including those 
with only one star) in our fit.'^ The resulting parame- 
ter values are displayed in Table V. Here the resonance 
masses are fixed at the respective (centroid) values given 
in Ref. [20]. The resulting resonance widths are all con- 
sistent with those quoted in Ref. [20]. The Pii(1440) 
resonance has practically no influence on the observables 
considered here and, therefore, it has been omitted in 



^ There are also established spin-5/2, -7/2, and -9/2 resonances 
exception of the latter resonance, all the resulting res- [20] in the energy region covered by the JLab data, but they 

onance masses are consistent with those quoted by the have been omitted in the present analysis. 
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TABLE III: Same as 
here to see whether 
X'/N = 1.04. 



Table II. More resonances were added 
this would further improve the fit. 
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FIG. 5: (Color online) Same as in Fig. 3 for the fit result of 
Table III. 



the fit set shown. For the NN-q' coupling constant, we 
obtained gNN-q' = 1-9. 

All these parameter sets illustrate the fact that cross 
sections do not impose enough constraints to the fit in 
order to extract definitive information on the resonances. 
Spin observables, on the other hand, do impose more 
stringent constraints and help distinguish among these 
parameter sets, as we shall show later. 

Although the parameter sets in Tables I-V yield com- 
parable fits to the cross section, the corresponding dy- 
namical contents are quite different from each other. 
Let us discuss, therefore, some of the different features 
present in the results corresponding to the various param- 
eter sets. Figure 3 shows some details of the dynamical 
content of our model corresponding to the fit results given 
in Table I. Here, both the Sn and P13 resonances have 
the largest contribution at low energies; the former dies 
out as the energy increases while the latter contribution 
persists to higher energies. The angular shape of the P13 
resonance current contribution is concave with a maxi- 
mum at 0,y m 90°. The Pn resonance contributes mostly 
around — 1.879 GeV. Its angular shape is rather flat 
(note that it includes both the s- and u-channel contri- 
butions). However, its interference with other contribu- 
tions, such as that due to the resonance, leads to a 
distinctive angular dependence. Although the Di^ res- 
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TABLE IV: Same as Table I. No P13 resonance was allowed 
here. x^/N = 1.10. 
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onance current is relatively small in this particular fit 
set, it plays an important role in reproducing the data 
through its interference with other currents. The mesonic 
current contribution plays a crucial role in reproducing 
the observed forward-peaked angular distribution, espe- 
cially at higher energies. This is a general feature ob- 
served in many reactions at high energies where the t- 
channel mechanism (either Regge trajectories or meson 
exchanges) accounts for the small-t behavior of the cross 
section. However, the present result shows also a compet- 
ing mechanism due to resonances and that the observed 
forward-peaked angular distribution is a result of signifi- 
cant interference effects. We note that this feature is not 
restricted to the particular set of the parameter values 
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cos(9^,) cos(0^,) cos(9^,) cos(6^,) 



FIG. 6: (Color online) Same as in Fig. 3 for the fit result of 
Table IV. The nucleonic current contribution (long-dashed 
curves) is negligible and cannot be seen on the present scale. 



of Table I, but it is also found in other sets that fit the 
data (note, in particular, a relatively large resonance 
contribution in Fig. 5 and a D13 resonance contribution 
in Figs. 4, 6 and 7 at higher energies). Therefore, this 
feature prevents us from fixing uniquely the mesonic cur- 
rent from the cross-section data at forward angles and 
higher energies. The nucleonic current contribution is 
very small here; however, as mentioned in the beginning 
of this section, the cross-section data alone do not im- 
pose stringent constraints on the fit so that it is possible 
to reproduce the data equally well with a much larger 
coupling constant, as can be seen in Tables III and V. 
We will come back to this issue later, in Sec. IIIC. 

Figure 4 shows the dynamical content of our model 
corresponding to the fit results given in Table II. Here, 
the 5*11 resonance contribution is considerable only at 
the lowest energy measurement and already at = 
1.627 GeV, it becomes very small and it is practically 
negligible for higher energies. Both the Pu and P13 
resonance contributions exhibit similar features to those 
observed in Fig. 3, with both contributions reaching its 
maximum around T-y — 1.677 GeV. The D13 resonance 
contribution is large over the energy region considered, 
except at lower energies. Above T-y « 1.677 GeV, it is 
the largest contribution. Its angular shape changes dras- 
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TABLE V: Same as 
at their PDG values. 



Table I. All resonance masses are kept 
X'/N = 1.01. 
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FIG. 7: (Color online) Same as in Fig. 3 for the fit result of 
Table V. 



tically with energy, starting with a small negative curva- 
ture in the lower energy region and ending with a roughly 
convex shape in the higher energy region. Note that this 
energy dependence of the angular shape is due to an in- 
terference between the two D13 resonances with different 
masses. Although somewhat larger, the mesonic current 
contribution is essentially the same to that in Fig. 3. As 
has been pointed out above, there is a considerable inter- 
ference effects between the mesonic and the (D13) reso- 
nance currents at higher energies. The nucleonic current 
is practically zero in Fig. 4 since, as mentioned above, 
the resulting NNrj' coupling constant is very small. 

In the fit result of Table III shown in Fig. 5, the 
Sii resonance contribution is very strong especially in 
the lower energy region and is quite appreciable even at 
higher energies. The Pn resonance contribution basically 
shows the same feature as in the fit sets discussed above. 
The P13 resonance contribution exhibits a convex angular 
shape, just opposite to the concave shape shown in Figs. 3 
and 4. This difference is due to the relative sign difference 
between the coupling constants giNN^r]' and g2NN'ri' as 
compared to the results shown in Figs. 3 and 4. The I?i3 
resonance contribution is largest around — 1.829 GeV. 
The shape of the angular distribution is quite different 
from the other fit discussed above. Together with the 
P13 resonance, it describes some of the details of the ob- 
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served angular distribution around T-y = 1.779 GeV. The 
mesonic current is much larger in this fit than in the other 
fits. In particular, it largely overestimates the measured 
cross section in the higher energy region. Its destructive 
interference with the other currents brings down the to- 
tal contribution in agreement with the data. Once more, 
this shows that one has to be cautious in trying to fix the 
i-channel contribution using the cross-section data at for- 
ward angles and higher energies. The nucleonic current 
gives an appreciable contribution in this fit, especially 
at higher energies and backward angles due to the it- 
channel. The corresponding NNrj' coupling constant is 

gNNr)' = 1-33. 

Figure 6 shows the dynamical content of the fit re- 
sult of Table IV. The Sn resonance contribution is 
largest at the lowest energy, but it decreases quickly 
as the energy increases. The Pn resonance contribu- 
tion is largest around = 1.627 GeV, with more pro- 
nounced angular distribution than in the other fit results. 
The Di3 resonance contribution has a concave angular 
shape in the lower energy region and is largest at around 
= 1.829 GeV. For higher energies the angular shape 
changes and gives the largest contribution for forward an- 
gles apart from the mesonic current, the latter providing 
again the bulk of the observed rise of the cross section 
at forward angles. The P13 resonance is not included 
in this fit set. We note that, unlike in Fig. 5, some of 
the details of the observed angular distribution around 
T-y = 1.779 GeV is not well reproduced, indicating the 
importance of both the P13 and -D13 resonances. The 
nucleonic current contribution is practically zero. 

The dynamical content of the fit result given in Table V 
is shown in Fig. 7. Overall, except for the lower energies, 
the mesonic current yields the largest contribution. The 
Sii and P13 resonance contributions are important in 
the low energy region while the Pn, P13, and Di^ reso- 
nances are important in the higher energy region. The 
nucleonic current is non-negligible only for backward an- 
gles at higher energies. Here, one major difference from 
the other fit results is the rather pronounced bending 
downward of the cross section (solid curves) for forward 
angles at lower energies. Measurements of the cross sec- 
tions for more forward angles would tell us whether such 
a behavior would indeed be necessary. 



B. Total cross sections 

Figure 8 shows the predictions for the total cross sec- 
tions obtained by integrating the corresponding differen- 
tial cross sections shown in Figs. 3-7. Although these pre- 
dictions may suffer from considerable uncertainties due 
to differences in the corresponding differential cross sec- 
tions, especially at lower energies, they exhibit a common 
feature, i.e., the total cross sections (solid curves) seems 
to show a bump structure around W = 2.09 GeV which is 
caused mainly by the Pn and/or _Di3 resonance depend- 
ing on the fit set. Note that the PDG quotes a two-star 




1.9 2.0 2.1 2.2 
W (GeV) 

FIG. 8: (Color online) Total cross section for 7p — > prj' as 
a function of the total energy of the system, W = y/a. The 
panels from top to bottom correspond to the fit results of Ta- 
bles I-V, as indicated. The overall total cross sections (solid 
lines) are broken down according to their dynamical contri- 
butions, with line styles defined in Fig. 3. The two dashed 
vertical lines are placed to guide the eye through the two 
bump positions in all panels. 



£'13(2080) and a one-star Pii(2100) resonance which are 
practically at this bump position. There is also a one- 
star resonance, S'ii(2090), that is just at the bump and, 
therefore, might have contributed to its structure. How- 
ever, the angular distribution does not favor this possibil- 
ity. The total cross section also seems to exhibit a bump 
structure at a lower energy of around W = 1.96 GeV due 
to the Sii, Pii and/or P13 resonance depending on the 
fit set. The latter two resonances can also contribute to 
the broadening of this bump depending on the fit set, as 
can be seen in Fig. 8. A rather sharp rise of the cross sec- 
tion from the threshold is caused by the 6*11 resonance, 
except in the top two panels, where the P13 resonance 
also contributes to this rise. The structures exhibited 
by the total cross section, in particular the bump around 
W = 2090 GeV, are unlikely to be artifacts of the present 



11 



predictions and, consequently, we would expect them to 
show up in the actual total cross-section data when they 
are measured. 



C. NNt]' coupling constant 

As we have seen in this section, unfortunately the 
present analysis cannot determine the NNrj' coupling 
constant, since the available cross-section data can be 
reproduced equally well with different sets of parame- 
ters in which this coupling constant varies considerably. 
However, an upper limit of its value can still be esti- 
mated. One of the reasons why gNN-q' cannot be ex- 
tracted uniquely from the cross-section data is that the 
resonance currents, especially the one due to the Di^ res- 
onance, can give rise to the observed enhancement of the 
backward-angle cross section as shown in Figs. 4 and 6. 
Also, the Pii resonance current alone can lead to a fea- 
ture of the cross section similar to that due to the nucle- 
onic current, i.e., the enhancement of the backward-angle 
cross section at higher energies through the u-channel 
contribution. The resonance currents can also interfere 
destructively with the nucleonic current in which case 
one obtains a larger NNrj' coupling constant. In fact, in 
a very extreme case, we have obtained a fit value as large 
as gNNri' = 3.0. It is obvious, therefore, that a more un- 
ambiguous extraction of this coupling constant requires 
going to an energy region where the resonance contribu- 
tions are small. Figure 9 illustrates this point; here we 
show the fit result considering the data with energies at 
= 2.129 GeV and above only and assuming a scenario 
in which no resonance currents contribute at these ener- 
gies. The resulting fit parameters are gNNri' = 2.10 for 
the NNrj' coupling constant and A^, = 1264 MeV for the 
cutoff parameter in the form factor at the rj'v^ vertex. 
In Fig. 9 we see that the nucleonic and mesonic currents 
interfere with each other. However, the interference pat- 
tern is such that it does not cause any problem in fixing 
both the nucleonic and mesonic current contributions to 
a large extent. In any case, judging from the overall re- 
sults of our analysis, we would not expect gNNri' to be 
much larger than 2. 
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FIG. 9: (Color online) Fit result with no resonances. Qmnt)' 
2.10 and A„ = 1264 MeV. See caption of Fig. 4. 



D. Meson exchanges versus Regge trajectory 

It is well known that f-channel processes at high ener- 
gies (above ~3-5 GeV) may also be described by Regge 
trajectories. However, how far down in energy one can 
go with this description before an explicit inclusion of or- 
dinary meson exchange is required is still an open issue. 
We would expect a smooth transition from a description 
in terms of meson exchanges to one in terms of Regge 
trajectory as one goes up higher in energy. The issue of 
meson exchanges versus Regge trajectory is particularly 
relevant in the present context, for the extracted reso- 
nance parameters can depend on these two alternatives 
for modeling the i-channel contribution [10]. 

In their analysis of the SAPHIR data [12], Chiang et 
al. [21] advocate the use of Regge trajectories while other 
authors [22] have employed vector-meson exchanges. In 
our previous analysis of these data, we found [10] that 
they can be reproduced equally well using either meson 
exchanges (with form factors at the vq'^ vertices) or a 
Regge trajectory for the ^-channel contribution. It is in- 
teresting to see whether the same is true for the new 
CLAS data [11]. We have repeated the calculation with 
the Regge trajectory following Ref. [10]. In particular, 
we replace the t-channel meson exchange propagators by 
the corresponding Regge trajectories keeping everything 
else unchanged, except that the form factor at the vq'^ 
vertex is set to unity. All the free parameters of the 
model are refitted again. We found that the fit quality 
using the Regge trajectory is, at best, comparable to that 
obtained using the ordinary meson exchanges for the t- 
channel. For example, the corresponding to the fit set 
of Table II is x^/N = 1.17 compared to xV^ = 1-04 
obtained with explicit meson exchanges; for the fit set of 
Table V, the Regge-trajectory result is x^/N = 1.08, as 
compared to x^ jN = 1.01 using meson exchanges. Sim- 
ilar results are obtained for other fit sets considered in 
subsection HI. A. 



E. Spin observables 

We now turn our attention to spin observables. As 
we have shown in Fig. 2, cross sections do not impose 
very severe constraints on the model parameter values. 
We expect spin observables to be more sensitive in this 
respect. The predictions for the beam and target asym- 
metries corresponding to the fit results of Tables I-V are 
shown in Fig. 10. As we can see, unlike the cross sec- 
tions (see Fig. 2), the predictions vary considerably be- 
tween the different parameter sets. For energies where 
the beam asymmetry is less sensitive to the parameter 
sets, the target asymmetry is quite sensitive and vice 
versa. Therefore, overall, a combined analysis of these 
spin observables will impose much more stringent con- 
straints on the fit and should help determine better the 
model parameters. 
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cautious with such an identification of the resonances. 
As we have seen, the cross-section data alone do not im- 
pose enough constraints for an unambiguous determina- 
tion of the resonance parameters. In this connection, we 
have shown that the beam and target asymmetries can 
help impose more stringent constraints. Furthermore, 
there is a possibility that some of the resonances in the 
present work are mocking up background contributions, 
especially those due to the final-state interaction, which 
is not taken into account explicitly in our calculation. 
Obviously, effects of the final-state interaction should be 
investigated in future work before a conclusive identifi- 
cation of the resonances can be made. 

We have predicted a bump structure in the total cross 
section &t W k, 2.09 GeV (see Fig. 8). If this is con- 
firmed, the Di3(2080) and/or Pii(2100) resonance may 
be responsible for this bump. 

Our study also shows that the nucleonic current should 
be relatively small. However, contrary to the expectation 
in our earlier work [10], the new high-precision cross- 
section data do not allow to pin down this current contri- 
bution due to the possible presence of resonance currents, 
especially of the D13 resonance, which can also lead to an 
enhancement of the cross section for backward angles at 
higher energies, a feature that otherwise arises from the 
w-channel nucleonic current contribution. These compli- 
cations notwithstanding, assuming that for the very high 
end of the present data set resonance contributions can 
be neglected, we argue in Sec. Ill C that the upper limit 
of gNNri' can now be lowered to a value of gNNri' ^ 2, 
whereas our previous analysis [10] had suggested an up- 
per limit of gNNri' ^ 3. Further corroboration of this 
finding is needed. 

In this respect, it should be noted that the result per- 
taining here to the NNrf coupling constant is, of course, 
a model-dependent one. Indeed, what is relevant in our 
calculations is the product of gNN-q' and the associated 
hadronic form factor which accounts for the off-shellness 
of the intermediate nucleon. Moreover, our NNrj' cou- 
pling constant is defined at the on-mass-shcU point, i.e.. 
We have analyzed the new CLAS [11] data of the gNNrj' = Sat att;' (9^ = w^') while the coupling required in 
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FIG. 10: (Color online) Photon beam and target nucleon 
asymmetries Es (top panel) and Et (bottom), respectively, 
for 7p prj' as a function of the r/' emission angle 8^' in the 
center-of-momentum system. See the caption of Fig. 2 for the 
meaning of the diflferent curves. 



IV. SUMMARY 



7P — > prj' reaction within an approach based on a rel- 
ativistic meson-exchange model of hadronic interactions. 
The present model is an extension of the one reported in 
Ref. [10] and it includes the nucleonic and the mesonic 
as well as the nucleon-resonance currents. The latter in- 
cludes both spin-1/2 and -3/2 resonance contributions 
in contrast to our previous work [10], where only spin- 
1/2 resonances were considered. In addition, we employ 
energy-dependent resonance widths in the present work. 
The resulting reaction amplitude is fully gauge invariant. 

We have shown that the mesonic as well as the spin-1/2 
and -3/2 resonance currents are important to describe the 
existing data quantitatively. The observed angular dis- 
tribution is due to delicate interference effects between 
the different currents. In our analysis, most of the re- 
sulting resonances may be identified with known reso- 
nances [20]. We emphasize, however, that one should be 



Eq. (1) in connection with the origin of the nucleon spin is 
at g2 ^ 0. Since the rj' meson is a relatively heavy meson 
{rrijji w 957 MeV), we would expect that gNNri' [o^ = 0) 
will be considerably smaller than its value at = m^, 
because of the presence of the form factor which usually 
cut down the coupling strength. Therefore, we might well 
expect that the NNrj' coupling aX = {) to be negligibly 
small, consistent with zero. 

We have also shown that the mesonic current contri- 
bution cannot be fixed unambiguously from the existing 
cross-section data because of the possible presence of the 
resonance currents, especially the Z)i3 resonance. A pos- 
sibility to determine the i-channel current is to measure 
the cross sections at higher energies where the resonance 
contributions becomes negligible. 

Furthermore, we have found that using a Regge trajec- 
tory in the i-channel instead of explicit meson exchanges 
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yields overall fit qualities that arc, at best, comparable 
to those obtained with meson exchanges. This indicates 
that explicit p and u) exchanges, as employed here, are 
completely adequate to describe the t-channel degrees of 
freedom at the present energies. 

Finally, the results of the present work should provide 
useful information for further investigations, both exper- 
imentally and theoretically, of the 7A'^ — > Nrj' reaction. 
In particular, measurements of cross sections at smaller 
forward and larger backward angles than are available 
in the present data set would already help constrain the 
model parameters considerably, as can be seen in Fig. 2. 
Total cross sections should also be measured in order 
to confirm or dismiss the bump structures, especially 
aroimd W = 2.09 GeV, predicted in the present calcu- 
lation. In addition, it is expected that measurements of 
spin observables — such as beam and target asymme- 
tries shown in Fig. 10 — would impose more stringent 
limits on the range of permissible parameters and this 
would undoubtedly provide a much improved description 
of the resonances and their properties in the energy re- 
gion covered by the existing data. From the theoretical 
side, it is possible that the nuclcon resonances introduced 
in the present work are mocking up the background con- 
tributions not taken into account in the calculation. In 
this connection, it is extremely interesting to investigate 
effects of the final-state interaction which has not been 
treated explicitly in the present calculation. Unfortu- 
nately, at present no realistic model is available that can 
provide the relevant ri'N final-state interaction. In ad- 
dition, effects of higher-spin resonances that have been 
ignored in the present analysis should be investigated in 
the future. 
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APPENDIX A: SPIN-3/2 RESONANCE 
PROPAGATOR 

We employ here the Rarita-Schwinger (RS) choice for 
the free Lagrangian of a spin-3/2 particle with mass m, 

jC = Vi^A'^^V. , (Al) 

where 

A''- = -l[a''^(;^-m)]+ , (A2) 

with p = id, the anticommutator bracket [a, 6]+ = ab + 
ba, and af" = ^{YY - Yl^)- (In the RS choice, the 



parameter A that usually appears in C is taken as A = —1 
[23].) From 

^i^xS^" = S^xA^- = g; , (A3) 
the propagator is then found as 

= »±ir^)^ = ^::;a±ji , ,a4) 

p-^ — p-^ — 

where A'''' is the RS tensor of Eq. (7) (with m = uir) 
and 

A^-^ = -gf^- + + ^ - ^ ^ g J ^ , (A5) 

which differs from (7) by the sign of the last term. 

When seeking an ansatz for describing a spin-3/2 res- 
onance, wc note first that there arc, of course, infinitely 
many ways to achieve a pole description whose on-shell 
behavior on the real axis corresponds to replacing the 
mass of the elementary propagator by 

m —> niR - i^ , (A6) 

where m.R is the resonance mass and Tr the associated 
width. In constructing a resonant propagator, we are 
guided by the following motivation. As in the spin-1/2 
case of (5), we want to describe the spin in terms of 
the elementary operators, i.e., we want to preserve the 
numerator structure of (A4) and the symmetry between 
the RS tensors A and A. In a schematic matrix notation, 
we therefore make the ansatz 

S = X-\^ + m)A = A{^ + m)X-'^ , (A7) 

putting, in analogy to the denominator of the spin-1/2 

case (5), 

X = - m|j)5 + iAij) + tur) , (A8a) 
X = (p2 - m\)g + iij) + mR)A , (A8b) 

with the operators A and A to be determined such that 
the second equality in (A7) holds true, i.e., 

X-^'iJ) + m)A = + m)X-^ . (A9) 

Multiplying this equation by A from both sides, one im- 
mediately finds the condition 

ii, + mR)Ak = KA{i) + tur) . (AlO) 

In view of Eq. (A3) and the fact that on-shell, at p^ = m\ 
and acting on a spin-3/2 eigenstate, the propagator must 
provide the width information, we find that the ansatz 

A^"" = -A^"'^ and A^"" = -A'"'^ , (All) 
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satisfies all constraints. T here may be any conveniently 
chosen width function that goes to the static width 
at the resonance mass mR. We thus have 



S{p) = 



- mR)g - i^T 



or 



Sip) = A 



A 

- mR)g - i—T 



(A12a) 



(A12b) 



By construction, both forms are completely equivalent, 
similar to the equivalence of both forms for the elemen- 
tary propagator (A4). 

The inversion here is to be performed on the full 16- 
dimensional space of Lorentz indices and component in- 
dices. There are various equivalent ways to do this; we 
have done it by introducing indices 



4/i + (3 and j — + a 



(A13) 



where /x, = 0, 1, 2, 3 are the Lorentz indices and (3,a = 
1, 2, 3, 4 are the component indices, and defining 16 x 16 



numerator and denominator matrices by 



+ ^R 



and 



Dij = Pxjpa^n,, - mRSn^6f3a - iNij- , (A15) 

respectively. Numerically inverting the denominator ma- 
trix D, we then calculate the spin-3/2 propagator as 



(A16) 



where k = Ap + a and summation over p is implied, as 
usual. 
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